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SUMMARY 


This  final  report  for  Contract  F33615-77-C-5093  summarizes  the  work  performed  during  the 
period  1  September  1977  to  30  June  1979. 

The  fatigue  crack  growth  behavior  of  an  advanced  powder  metallurgy  (P/M)  superalloy 
(AF2-1DA)  was  characterized,  and  an  interpolative  mathematical  model  was  put  into  computer 
code  to  permit  accurate  predictions  of  crack  propagation  rates  under  different  stress- 
temperature-time  conditions.  The  Interpolative  Hyperbolic  Sine  Model  (References  1,  2  and  3) 
having  the  general  equation 

log  (da/dN)  =  C,  sinh  (C2(log  (AK)  +  Cs))  +  C4  (1) 

was  the  basis  for  model  development. 

Interpolative  crack  growth  rate  models  are  presented  showing  the  effects  of  frequency,  stress 
ratio,  temperature,  and  dwell  time  at  maximum  tensile  stress.  The  effect  of  specimen  thickness 
on  crack  growth  rates  is  also  shown. 

Verification  testing  to  check  the  model  demonstrated  the  model’s  capability  of  predicting 
crack  growth  rates  where  no  data  exist.  Crack  growth  rates,  calculated  using  the  interpolative 
model  before  testing,  were  used  to  predict  specimen  cyclic  lives.  The  results  of  these 
predictions  (Figures  70  through  75)  show  that  the  interpolative  hyperbolic  sine  (SINH)  model 
accurately  predicts  crack  propagation  life  of  AF2-1DA  at  elevated  temperatures. 
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SECTION  I 
INTRODUCTION 


The  presence  of  a  crack  in  a  stressed  component  necessitates  redistribution  of  stresses 
around  the  crack.  The  stress  intensity  factor  is  a  parameter  that  reflects  this  redistribution  and 
is  a  function  of  nominal  stress,  flaw  size,  and  specimen  and  crack  geometries.  The  concept  of 
stress  intensity  factor  was  originally  defined  for  an  infinitely  sharp  crack  in  a  perfectly  elastic 
medium.  In  most  engineering  materials,  localized  plastic  deformation  occurs  due  to  high  stresses 
at  the  crack  tip,  and  it  is  this  deformation  that  gives  the  materia)  resistance  to  crack  propagation. 
In  a  completely  brittle  material,  crack  tip  stress  relaxation  is  so  small  that  simple  reinitiation  of 
a  stopped  crack  is  sufficient  to  promote  complete  fracture. 

The  degree  of  brittleness  of  a  material  (and  the  limit  to  the  applicability  of  linear  elastic 
fracture  mechanics)  is  directly  related  to  the  type  of  relaxation  process  that  occurs  at  the  crack 
tip.  In  the  high  temperature  fatigue  process,  this  relaxation  is  expected  to  depend  on  the  relative 
degree  of  elastic,  plastic,  creep,  and  chemical  work  expended  at  the  crack  tip.  In  this  program, 
linear  elastic  fracture  mechanics  is  investigated  only  to  the  extent  required  to  determine 
specimen  thickness  needed  for  thickness-independent  data.  This  is  adequate  because  only  the 
compact  specimen  is  being  used  and  anomalies  associated  with  high  net  section  stresses  and 
geometric  variations  are  avoided. 

The  technical  approach  of  this  program  consists  of  (1)  isolation,  testing  and  modeling  the 
effects  of  stress,  time  and  temperature  on  the  crack  growth  rate  of  the  advanced  P/M  superalloy 
AF2-1DA;  (2)  synthesis  and  testing  of  models  that  accurately  reflect  parametric  interactions;  and 
(3)  demonstration  of  the  models  for  a  simple  crack  geometry  subjected  to  a  typical  stress-time- 
temperature  history. 

The  technical  program  consists  of  four  phases:  The  first  phase  determines  the  effect  of 
specimen  thickness  on  crack  growth  rates  for  AF2-1DA  material  and  characterizes  the  effects  of 
stress  ratio,  cyclic  frequency,  and  temperature.  Phase  II  consists  of  generalized  interpolative 
model  development  based  on  the  hyperbolic  sine  equation 

log  (da/dN)  =  C,  sinh  (C2  (log  (AK)  +  C,))  +  C«  (1) 

and  Phase  HI  provides  for  demonstration  of  the  predictive  model.  Phase  IV  encompasses  the 
reporting  requirements. 
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TECHNICAL  PROGRESS 


MATERIAL  PROCUREMENT 

This  program  examined  the  crack  growth  behavior  of  an  advanced  powder  metallurgy 
(P/M)  superalloy,  AF2-1DA.  This  alloy  was  developed  as  a  turbine  disk  material  for  use  at 
1400°F,  a  temperature  at  which  time  dependent  deformation  is  a  major  concern.  The  material 
was  processed  by  the  GATORIZING®  technique  of  superplastic  forming,  creating  the  micro¬ 
structure  of  a  typical  disk  forging.  All  testing  was  performed  on  a  single  heat  of  material.  Low- 
cycle  fatigue  (LCF),  creep,  tensile,  and  combination  stress  rupture  tests  were  performed  to 
qualify  the  material.  The  specimens  are  shown  in  Figures  1  through  4. 

The  first  forging  developed  a  maze  of  cracks  during  heat  treatment  to  AMS  5881.  It  was 
determined  that  this  specification  contained  an  error  in  the  heat  treatment  procedure.  The 
solution  heat  treatment  calls  for  an  oil  quench  and  it  should  require  rapid  air  cooling.  The 
material  used  in  this  program  was  heat-treated  by  taking  exception  to  AMS  5881  and  the 
microstructure  compared  favorably  to  a  good  full-scale  component  forging,  an  F100  3rd-stage 
turbine  disk  (S/N  455-B2)  developed  under  AFML  Contract  F33615-74-C-5108.  Figure  5  shows 
the  structures  to  be  identical  in  all  important  respects.  Comparison  of  the  size  and  distribution 
of  y'  reveals  that  the  heat  treatment  developed  under  AFML  Contract  F33615-74-C-5108 
(reference  4)  was  duplicated.  The  grain  sizes  of  the  components  were  also  similar,  ranging  from 
ASTM  0  to  3,  predominantly  0  and  1. 

The  qualification  testing  for  the  three  forgings  is  summarized  in  Table  1.  AMS  5881  requires 
the  stress  rupture  specimen  to  rupture  in  not  less  than  23.0  hr,  but  the  test  specimen  lasted 
between  15.7  and  19.3  hr  (except  for  the  first  BAUD-3  specimen  that  failed  early  in  the  notch). 
Out  of  six  disks  tested  for  AFML  Contract  F33615-74-C-5108,  only  one  would  have  met  the  AMS 
5881  stress  rupture  requirements  for  time  to  rupture.  Data  from  the  six  disks  show  an  average 
stress-rupture  life  of  18.9  hr  and  a  97.5%  lower  bound  of  11.5  hr  based  on  12  data  points  and  an 
assumed  log-normal  distribution  of  life.  Tests  performed  on  the  AF2-1DA  forgings  resulted  in 
stress-rupture  lives  near  the  18.9  hr  mean  life.  This  indicates  material  behavior  similar  to  the 
results  reported  in  the  referenced  contract. 

The  tensile  properties  recommended  in  AFML-TR-76-101  are  also  reduced  from  AMS  5881. 
Table  I  summarizes  the  tensile  requirements  of  both  specifications  in  addition  to  the  test  results. 

Four  (4)  isothermal,  axial  strain-controlled  LCF  tests  were  performed  to  provide  a 
comparison  of  the  LCF  capability  of  the  forging  used  for  test  material  in  this  program  with 
previously  tested  AF2-1DA.  All  tests  were  conducted  at  760°C  (1400°F)  at  a  cyclic  frequency  of 
0.17  Hz  (10  cpm)  using  a  sawtooth  strain  vs  time  waveform.  The  strain  cycle  was  all  tensile  such 
that  the  mean  strain  was  equal  to  one-half  the  maximum  strain  (R<  =  0). 

The  AF2-1DA  material  used  in  this  program  exhibited  LCF  lives  comparable  to  previously 
tested  AF2-1DA.  The  test  results  are  presented  in  Table  2  and  plotted  in  Figure  6  for  comparison. 
The  total  strain  range,  elastic  and  inelastic  strain  components,  stress  range,  and  mean  stress  were 
determined  at  specimen  half-life  (Nr/2)  and  are  presented  in  Table  2.  All  four  specimens 
exhibited  stable  cyclic  stress-strain  behavior.  Failure  origins  for  specimens  tested  at  the  higher 
total  strain  ranges  (1.1%  and  0.9%)  were  surface  connected;  the  failure  origin  in  the  specimen 
tested  at  0.82%  strain  range  was  subsurface  and  at  a  probable  nonmetallic  inclusion.  The  test 
specimen  cycled  at  0.6%  total  strain  range  was  discontinued  after  182,000  cycles  with  no 
indications  of  failure. 
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Figure  2.  Creep  Specimen 
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NOTES: 

1.  Ml  dimensions  In  inches. 
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to  within  .0005. 

6.  Identification  markings  permitted  only  on  specimen  ends. 

7.  per  Peterson,  “Stress  Concentration  Factors". 


Mag:  10.000X 

B.  Microstructure  of  F100  3rO*Stage  Turbina  Disk  Devalopad 
Under  AFML  Contract  F3361 5-74-C-5108 


Figure  5.  Comparison  of  AF2-1DA  Pancake  Microstructure  to 
the  Microstructure  of  FlOO  3rd-Stage  Turbine  Disk 
from  AFML  Contract  F33615-74-C-5108 


TABLE  1. 

QUALIFICATION  TEST  RESULTS 


Tensile  Properties 


Tensile 


Temp. 

CC) 

0.2 %  Y.S. 
(MPa) 

Strength 

(MPa) 

El 

% 

Ra 

% 

AMS  5881t 

R.T. 

965 

1310 

10.0 

12.0 

AFML-TR-76-101t 

R.T. 

931 

1310 

10.0 

12.0 

BAUD-3 

R.T. 

981 

1517 

14.5 

15.6 

BAUD-4 

R.T. 

960 

1480 

14.5 

16.3 

BAUD-5 

R.T. 

983 

1488 

16.0 

17.5 

AMS  5881f 

816 

862 

1034 

12.0 

AFML-TR-76-lOlf 

816 

1000 

827 

■m 

BAUD-3 

816 

861 

985 

28.7 

BAUD-4 

816 

907 

1056 

18.0 

29.0 

BAUD-5 

816 

848 

994 

24.4 

Creep  Properties 


Heat 

Code 

Temp. 

(°C) 

Stress 

(MPa) 

Rupture 

Time 

(Hours) 

Hours  To  Creep 
0.1%  0.2% 

%  at 
100  Hrs 

BAUD-3 

760 

483 

>100 

12.0 

55.6 

.312 

BAUD-4 

760 

483 

>100 

53.0 

68.7 

.273 

BAUD-5 

760 

483 

>100 

31.6 

90.6 

.215 

Stress 

Rupture 

Properties 

Rupture 

Heat 

Temp. 

Stress 

Time 

El 

RA 

Code 

CC) 

(MPa) 

(Hours) 

% 

% 

BAUD-3 

816 

552 

9.0* 

_ 

_ 

BAUD-3 

816 

552 

15.7 

8.1 

12.4 

BAUD-3 

816 

552 

19.3 

7.7 

11.6 

BAUD-4 

816 

552 

17.5 

7.6 

10.8 

BAUD-5 

816 

552 

16.5 

6.9 

8.0 

+ Requirements  of  this  specification 

*  Failed  in  notch  after  9.0  hr,  continued  to  test  smooth  section 
with  failure  after  19.4  hr. 


TABLE  2. 


AXIAL  CONTROLLED  STRAIN  LOW-CYCLE  FATIGUE  RESULTS  OF 
AF2-1DA  TESTED  IN  AIR  760°C  (1400°F),  0.17  Hz  (10  cpm),  R,  =  0. 


Cycles 

To 


Spec. 

S/N 

Strain  (m/m  at  N,/2) 

Mean  Stress  . 

Stress  Range 

Cyclic 

Stability 

Failure 

Range 

Elastic 

Inelastic 

Creep 

at  St/2 

Cycle  I 

N,/2 

N, 

1' 

0.0112 

0.0106 

0  0006 

0 

119  MPa 
(17.3  kai) 

1737  MPa 
(252.0  ksi) 

1726  MPa 
(250.4  ksi) 

Stable 

1,037 

2' 

0.0090 

0.0088 

0.0002 

0 

148  MPa 
(21.4  ksi) 

1476  MPa 
(214.1  ksi) 

1442  MPa 
(209.2  ksi) 

Stable 

5,282 

3" 

0.0082 

0.0081 

0.0001 

0 

267  MPa 
(38.7  ksi) 

1306  MPa 
(189.4  ksi) 

1317  MPa 
(191.0  ksi) 

Stable 

6,623 

4". 

0.0060 

0.0060 

<0.0001 

0 

108  MPa 
(15.7  ksi) 

988  MPa 
(143.3  ksi) 

966  MPa 
(140.1  ksi) 

Stable 

182,460 

'Surface  origin. 

"Origin  approximately  0.4  mm  from  specimen  surface  at  possible  inclusion. 
"'Test  discontinued.  Did  not  fail. 


10*  103  104  10® 

Cycles  to  Failure  ro  issue 

Figure  6.  Axial  Controlled  Strain  Low-Cycle  Fatigue  Crack  Initiation  Results  for 
AF2-1DA  Tested  in  Air  at  760°C  (1400°F),  0.17  Hz  (10. cpm),  R,  *  0 


Figure  7  compares  the  results  of  a  12.7  mm  (0.500  in.)  thick  compact  specimen  tested  at 
650°C  (1200°F),  R  =  0.1, 10  cpm  with  data  from  previous  experience.  The  excellent  agreement  of 
the  data  with  previous  experience  removes  all  doubt  as  to  whether  this  heat  of  material  is  typical 
of  other  AF2-1DA  material. 

TEST  SPECIMEN  AND  PROCEDURES 

Figure  8  shows  the  compact  specimen  used  to  obtain  crack  propagation  data.  Testing  was 
conducted  on  servohydraulic,  closed-loop,  load-controlled  testing  machines.  Specimens  were 
precracked  using  procedures  outlined  in  ASTM  E-399.  Crack  lengths  were  measured  directly  with 
a  Gaertner  traveling  microscope  without  cooling  the  specimen.  The  intervals  between  crack 
growth  measurements  were  selected  to  obtain  crack  growth  increments  of  approximately  0.5  mm 
(0.02  in.),  which  normally  results  in  an  average  of  40  to  50  readings  per  specimen.  The  crack 
propagation  tests  were  conducted  with  a  triangular  loading  wave  form  or  a  triangle  wave  with  a 
hold  time  at  the  tensile  peak  with  all  portions  of  the  cycle  under  tensile  load -controlled 
conditions.  Table  3  lists  all  crack  propagation  test  specimens  for  AF2-1DA. 
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da/dN,  in. /Cycle 


da/dN,  mm/Cycle 


TABLE  3. 

CRACK  PROPAGATION  TEST  SPECIMENS  FOR  AF2-1DA 


Specimen 

Number 

Temperature 

(°C)  CF) 

Cyclic 

Frequency 

Stress 

Ratio 

Thickness 

(in.) 

Remarks 

1101 

649 

1200 

0.17  Hz  (10  cpm) 

0.10 

0.500 

1102 

427 

800 

0.17  Hz  (10  cpm) 

0.10 

0.503 

1103 

649 

1200 

120  Sec.  Dwell 

0.10 

0.502 

1104 

649 

1200 

0.0063  Hz  (0.5  cpm) 

0.10 

0.500 

1106 

649 

1200 

0.17  Hz  (10  cpm) 

0.50 

0.468 

1107 

427 

800 

0.0063  Hz  (0.5  cpm) 

0.10 

0.500 

Nc  data 

1108 

649 

1200 

120  Sec.  Dwell 

0.10 

0.500 

No  data 

1109 

427 

800 

20  Hz 

0.10 

0.487 

1110 

760 

1400 

0.0083  Hz  (0.5  cpm) 

0.10 

0.501 

Overstress  failure 

1111 

760 

1400 

120  Sec.  Dwell 

0.10 

0.501 

1112 

427 

800 

0.17  Hz  (10  cpm) 

0.50 

0.250 

1114 

760 

1400 

600  Sec.  Dwell 

0.10 

0.498 

Excessive  Curvature 

1115 

649 

1200 

0.0083  Hz  (0.5  cpm) 

0.50 

0.500 

1117 

760 

1400 

120  Sec.  Dwell 

0.10 

0.124 

Equipment 

1118 

649 

1200 

0.17  Hz  (10  cpm) 

0.10 

0.125 

1119 

427 

800 

0.17  Hz  (10  cpm) 

0.10 

0.124 

Equipment  Failure 

1120 

427 

800 

0.17  Hz  (10  cpm) 

0.10 

0.250 

1122 

760 

1400 

120  Sec.  Dwell 

0.10 

0.124 

1124 

649 

1200 

600  Sec.  Dwell 

0.10 

0.452 

1125 

649 

1200 

20  Hz 

0.10 

0.249 

1126 

649 

1200 

0.17  Hz  (10  cpm) 

0.10 

0.249 

1128 

649 

1200 

120  Sec.  Dwell 

0.10 

0.130 

1129 

760 

1400 

0.17  Hz  (10  cpm) 

0.10 

0.124 

1130 

649 

1200 

120  Sec.  Dwell 

0.10 

0.254 

1131 

427 

800 

0.17  Hz  (10  cpm) 

0.10 

0.126 

1132 

427 

800 

20  Hz 

0.10 

0.253 

Equipment  Failure 

1134 

427 

800 

20  Hz 

0.10 

0.243 

1135 

649 

1200 

20  Hz 

0.50 

0.244 

1136 

760 

1400 

20  Hz 

0.50 

0.244 

1137 

427 

800 

0.0063  Hz  (0.5  cpm) 

0.50 

0.244 

1138 

427 

800 

20  Hz 

0.10 

0.257 

1139 

427 

800 

20  Hz 

0.80 

0.254 

1140 

649 

1200 

0:17  Hz  (10  cpm) 

0.80 

0.250 

1142 

649 

1200 

20  Hz 

0.80 

0.249 

1143 

760 

1400 

0.0083  Hz  (0.5  cpm) 

0.50 

0.248 

1144 

649 

1200 

0.0083  Hz  (0.5  cpm) 

0.80 

0.249 

1145 

760 

1400 

120  Sec.  Dwell 

0.10 

0.250 

1146 

427 

800 

20  Hz 

0.50 

0.250 

1147 

760 

1400 

0.0083  Hz  (0.5  cpm) 

0.80 

0.249 

1148 

760 

1400 

120  Sec.  Dwell 

0.10 

0.500 

1149 

760 

1400 

0.17  Hz  (10  cpm) 

0.10 

0.501 

1150 

649 

1200 

Sustained  Load 

1.00 

0.500 

1151 

760 

1400 

Sustained  Load 

1.00 

0.498 

1152 

649 

1200 

300  Sec.  Dwell 

0.10 

0.499 

1153 

718 

1325 

15  Sec.  Dwell 

0.10 

0.499 

Verification  Test 

1154 

718 

1325 

15  Sec.  Dwell 

0.10 

0.500 

Verification  Test 

1156 

760 

1400 

20  Hz 

0.10 

0.250 

1157 

760 

1400 

0.0083  Hz  (0.5  cpm) 

0.10 

0.252 

1160 

760 

1400 

0.17  Hz  (10  cpm) 

0.80 

0.251 

1162 

649 

1200 

0.17  Hz  (10  cpm) 

0.80 

0.250 

1165 

649 

1200 

30  Sec.  Dwell 

0.10 

0.300 

1167 

538 

1000 

0.017  Hz  (1  cpm) 

0.65 

0.242 

Verification  Test 

1168 

538 

1000 

0.017  Hz  (1  cpm) 

0.65 

0.251 

Verification  Test 

1173 

760 

1400 

600  Sec.  Dwell 

0.10 

0.750 

DATA  ANALYSIS  PROCEDURES 

The  reduction  of  crack  length,  a,  versus  cycles,  N,  data  is  accomplished  by  using  the 
incremental  polynomial  method  (Reference  5).  This  method  for  computing  da/dN  involves  fitting 
a  second-order  polynomial  (parabola)  to  sets  of  seven  successive  data  points.  The  equation  for  the 
local  fit  takes  the  form: 


=  b0  +  b. 


rN‘-H+b, 

r  n4-i, 

-V 

L  &  J 

-  £2 

j 

where: 


-1  < 


<  +  1 


(2) 


and  b„,  b,  and  b2  are  the  regression  parameters  which  are  determined  by  the  least  squares  method 
(that  is,  minimization  of  the  square  of  the  deviations  between  observed  and  fitted  values  of  crack 
length)  over  the  range  ai  to  a7.  The  value  4,  is  the  fitted  value  of  crack  length  at  N4.  The 
parameters,  £,  =  Vt  (N,  +  N,)  and  =  Vi  (N,  -  N,),  are  used  to  scale  the  input  data.  The  crack 
growth  rate  at  N,  is  obtained  directly  from  the  derivative  of  equation  2  as  follows: 

(da/dN)  d,  =  ^  +  2  b2  (N4  -  £,)  /  g  (3) 

The  AK  associated  with  this  da/dN  is  calculated  using  the  fitted  crack  length^,. 

After  calculation  of  da/dN  and  AK  using  the  first  seven  a  versus  N  data  points,  another 
da/dN  versus  AK  value  is  calculated  using  actual  a  versus  N  data  points  2  through  8  (renumbered 
1  through  7  for  use  in  the  equations).  The  local  fit  is  moved  along  the  data  one  point  at  a  time 
until  the  last  seven  points  are  used  in  the  calculation.  A  computer-drawn  plot  of  actual  crack 
length  versus  cycles  is  produced  when  the  da/dN  and  AK  calculations  are  made  to  allow  the 
actual  data  to  be  scrutinized  by  the  engineer. 

The  da/dN  versus  AK  calculations  from  the  seven-point  incremental  polynomial  computer 
program  are  modeled  using  the  hyperbolic  sine  equation,  and  the  final  da/dN  versus  AK 
computer-drawn  plots  are  produced  on  log-log  scale.  Families  of  da/dN  versus  AK  curves  are  then 
related  by  the  SINH  coefficients,  C,  through  C4,  to  provide  frequency,  stress  ratio,  and 
temperature  models. 

The  final  interpolative  models  are  presented  using  data  reduced  by  the  incremental 
polynomial  method;  however,  many  of  the  preliminary  models  presented  data  evaluated  by  the 
direct  secant  method.  Discrete  values  of  Aa  and  AN  were  computed  from  raw  laboratory  data.  By 
not  smoothing  (regressing)  the  a,  N  data,  the  actual  local  Aa/AN  perturbations  are  observable  in 
the  final  da/dN  versus  AK  curves.  Figures  9  through  15  compare  the  direct  secant  and 
incremental  polynomial  methods,  and  the  only  significant  difference  for  these  well-behaved  data 
is  three  additional  data  points  at  each  extreme  for  the  direct  secant  method. 

Crack  propagation  under  constant  amplitude  loading  conditions  is  a  function  of  the  applied 
stress  intensity  range  (within  the  limits  of  applicability  of  tinear  elastic  fracture  mechanics).  The 
applied  stress  intensity  range,  AK,  is  the  driving  force  for  crack  propagation.  Many  relationships 
have  been  developed  to  correlate  observed  crack  growth  rate  and  stress  intensity.  Paris  presented 
the  simple  relationship  (Reference  6): 


da/dN  =  C  (AK) 


(4) 


where  C  and  n  are  material  constants.  At  elevated  temperatures,  however,  the  crack  growth 
process  is  a  complicated  function  of  stress  ratio,  temperature,  load  history,  and  environment. 
These  dependencies  make  the  general  use  of  equations,  such  as  equation  (4),  more  difficult.  A 
model  developed  at  P&W A/Florida  (References  1,  2  and  3)  was  used  to  describe  the  effects  of 
cyclic  frequency,  stress  ratio,  and  temperature  on  the  crack  growth  rate  of  AF2-1DA.  The  model 
is  based  on  the  hyperbolic  sine  equation: 

log(da/dN)  =  C,  sinh  (Cs(log  (AK)  +  C,))  +  C4,  (1) 

where  the  coefficients  have  been  shown  to  be  functions  of  test  frequency,  stress  ratio,  and 
temperature: 

C-  =  material  constant 

Cs  =  f,  (R,  *,  T) 

C3  =  fa  (C41  V,  R) 

C4  =  f4  (v,  R,  T). 

The  hyperbolic  sine  equation  was  selected  as  the  model  for  the  following  reasons: 

•  It  exhibits  the  overall  shape  of  typical  da/dN  vs  AK  plots  obtained  over 
several  decades  of  crack  growth  rates. 

•  All  or  part  of  the  equation  may  be  used  to  fit  data  since  the  hyperbolic  sine 
has  both  a  concave  and  a  convex  half  and  a  nearly  linear  portion  near 
inflection.  Also,  the  slope  at  inflection  can  vary  with  the  fitting  constants.  (By 
comparison,  the  slope  of  an  x3  model  is  always  zero  at  inflection.) 

•  The  sinh  is  not  periodic  (e.g.,  trigonometric  tangent)  nor  asymptotic  (e.g., 
tangent,  or  inverse  hyperbolic  tangent);  therefore,  when  extrapolation 
becomes  necessary,  the  sinh  behaves  well  at  distances  removed  from  the  data, 
quite  unlike  most  polynomials,  periodic,  or  asymptotic  functions. 

•  This  model  requires  no  information  other  than  a,  N  data.  By  comparison, 
some  other  models  in  current  use  require  both  Kth  and  KJC.  in  addition  to  a. 

N  data,  to  model  crack  growth  behavior.  Both  K,h  and  K,c  are  difficult  to 
obtain  experimentally;  K,h  because  of  the  extremely  small  crack  growth 
measurements  necessary,  and  K,c  because  of  gross  plasticity  at  the  crack  tip 
encountered  in  fracture-toughness  testing  at  elevated  temperatures. 

The  hyperbolic  sine  is  defined  as: 

y  =  sinh  x  =  - — ~ —  (5) 

and  when  presented  on  Cartesian  coordinates,  it  appears  as  shown  in  Figure  16.  The  function  is 
zero  at  x  =  0  and  has  its  inflection  there. 

The  introduction  of  the  four  regression  coefficients,  Ci  through  C4,  permits  relocation  of  the 
point  of  inflection  and  scaling  of  both  axes.  In  the  equation, 

(y  ~  C4)  =  sinh  (x  +  C,),  (6) 

C,  establishes  the  horizontal  location  of  the  hyperbolic  sine  point  of  inflection  and  C4  locates  its 
vertical  position. 
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Figure  11.  Comparison  of  Direct  Secant  and  Incremental  Polynomial  Data 
Reduction  Methods  for  Specimen  1103 
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Figure  12.  Comparison  of  Direct  Secant  and  Incremental  Polynomial  Data 
Reduction  Methods  for  Specimen  1104 
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Figure  14.  Comparison  of  Direct  Secant  and  Incremental  Polynomial  Data 
Reduction  Methods  for  Specimen  1114 
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Figure  15.  Comparison  of  Direct  Secant  and  Incremental  Polynomial  Data 
Reduction  Methods  for  Specimen  1118 
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To  scale  the  axes,  C,  and  C2  are  introduced 


(y  -  C4) 
C, 


=  sinh  (C,  (x  +  C,)) 


(7) 


which  can  be  rewritten  as 

y  =  C,  sinh  (C2  (x  +  Ca))  +  C4  (8) 

of  which  equation  (2)  is  a  special  case  where  y  =  log(da/dN)  and  x  =  log(AK).  Note  that  C,  has 
units  of  log  (AK)  and  C4  has  units  of  log  (da/dN);  C,  and  C,  are  dimensionless  and  can  be 
conceptualized  as  stretching  the  curve  vertically  and  horizontally,  respectively.  Experience 
indicates  that  for  a  given  material,  can  be  fixed  without  adversely  affecting  model  flexibility 
(Reference  3). 

The  hyperbolic  sine  model  is  easily  adapted  to  describe  the  fundamental  parametric  effects 
of  stress  ratio,  frequency,  and  temperature  on  crack  growth  rate.  Experience  with  turbine  disk 
alloys  indicates  that  changing  test  frequency,  while  holding  stress  ratio  and  temperature 
constant,  produces  crack  growth  curves  similar  in  shape,  but  shifted  along  a  nearly  vertical  line 
passing  through  the  points  of  inflection.  The  location  of  these  inflection  points  is  related  to  test 
frequency.  Similar  results  can  be  obtained  with  stress  ratio  and  temperature.  Figure  17 
schematically  depicts  the  qualitative  effects  of  frequency,  stress  ratio,  and  temperature  on  crack 
growth  rates. 

Because  of  the  simple  relationships  observed  between  the  coefficients  of  the  SINH  model 
and  the  fundamental  propagation-controlling  parameters,  interpolations  are  straightforward.  It 
is  here  that  the  model  demonstrates  its  great  usefulness:  the  Hyperbolic  Sine  Model  provides 
descriptions  of  crack  propagation  characteristics  even  where  data  are  unavailable.  The 
interpolation  algorithm  is  described  in  Reference  7,  and  the  four  SINH  model  coefficients  can  be 
determined  as  follows: 

C,  =  material  constant  (0.5  for  AF2-1DA) 

C,  =  m2  log  (v)  +  b2 
Ca  =  m,  C4  +  ba 
C4  =  m4  log  (v)  +  b4. 

Computation  of  a  crack  growth  rate  equation  for  any  given  frequency  (with  stress  ratio  and 
temperature  held  constant)  is  a  straightforward  calculation,  once  the  above  linear  relationships 
have  been  established. 
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PhaM  I:  Crack  Growth  Data  Generation 

The  objectives  of  Phase  I  are  to  determine  the  effects  of  specimen  thickness  on  crack 
growth  rates  for  AF2-1DA  material  and  characterize  the  effects  of  stress  ratio,  cyclic  frequen¬ 
cy,  and  temperature. 

The  effects  of  thickness  on  crack  propagation  were  determined  by  testing  3.17  mm, 
6.35  mm,  and  12.70  mm  thick  specimens.  Thickness  effects  were  determined  at  427°C,  649°C, 
and  760°C  for  both  cyclic  and  cyclic/dwell  stress-time  waveforms. 

The  modeling  of  frequency  effects  on  subcritical  crack  growth  is  accomplished  using 
20  Hz,  0.17  Hz,  and  0.0083  Hz  as  the  test  points,  and  the  effects  of  stress  ratio  are  modeled 
using  data  generated  at  R  =  0.1,  R  =  0.5,  and  R=0.8.  The  temperature  effects  are  characterized 
at  isothermal  conditions,  and  the  modeling  of  these  effects  is  based  on  test  data  at  427°C, 
649°  C,  and  760°C. 

Applicability  of  Linear  Elastic  Fracture  Mechanics 

The  usefuless  of  linear  elastic  fracture  mechanics  (LEFM)  depends  on  a  uniparametrical 
relationship  between  crack  growth  rate  and  the  stress  intensity  factor.  Crack  tip  inelasticity, 
due  to  material  response  at  elevated  temperatures,  can  preclude  general  utility  of  K  as  the 
correlative  parameter.  Tests  were  conducted  to  ensure  sufficient  geometric  constraint  and/or 
environment  embrittlement  at  the  crack  tip  as  to  render  crack  growth  rate  —  AK  relationships 
independent  of  specimen  geometry  (e.g.,  only  thickness  in  this  case).  Universal  applicability  of 
LEFM  is  not  determined  in  this  program. 

In  advanced  P/M  alloys,  crack  tip  inelasticity  increases  with  increasing  applied  stress 
intensity.  When  the  plastic  zone  size  becomes  comparable  to  the  specimen  thickness,  the 
model  of  crack  propagation  changes  from  pure  Mode  I  (i.e.,  opening  mode)  to  mixed  mode 
(MODE  I  and  MODE  III,  tearing  mode)  and  LEFM  is  invalid.  For  a  constant  applied  stress 
intensity  range,  mixed  mode  propagation  rates  are  slower  than  those  for  plane  strain  constraint 
at  the  crack  tip.  Only  compact  specimens  were  used  in  this  program,  alleviating  problems  with 
differing  geometry  and  net  section  stress;  therefore,  the  only  requirement  was  that  a  specimen 
produce  thickness-independent  data  (i.e.,  a  thicker  specimen  will  not  produce  a  different 
propagation  rate). 

Figure  18  shows  the  effect  of  thickness  on  crack  growth  rates  for  3.18  mm  (0.125  in.), 
6.35  mm  (0.250  in.),  and  12.70  mm  (0.500  in.)  thick  specimens.  All  specimens  were  tested  at 
427°C  (800°F),  10  cpm,  R=0.1.  To  ensure  thickness-independent  crack  growth  data,  a 
minimum  specimen  thickness  of  6.35  mm  is  needed. 

The  cyclic  (0.17  Hz)  crack  growth  rates  at  649°C  (1200°F),  R  =  0.1  are  shown  in  Figure  19 
for  three  thicknesses  (3.18  mm,  6.35  mm  and  12.70  mm).  A  comparison  of  these  rates  shows 
that  6.35  mm  is  the  minimum  thickness  needed  for  thickness-independent  data  at  these 
conditions. 

At  760°C  (1400°F),  10  cpm,  and  R=0.1,  the  3.18  mm  thick  specimen  produced  nearly  the 
same  crack  growth  rate  (Figure  20)  as  the  12.70  mm  thick  specimen.  Since  the  two  extremes  in 
thickness  produce  nearly  equal  data,  6.35  ram  thick  specimens  are  probably  adequate  to 
produce  thickness-independent  data. 
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A  comparison  at  649°C,  R=0.1,  with  a  120  sec  dwell  at  maximum  tensile  load  (using  the 
same  ramp  as  a  0.17  Hz  cycle)  shows  that  even  the  12.70  mm  thick  specimen  may  not  produce 
thickness-independent  data.  However,  since  the  difference  between  the  6.35  mm  and  the 
12.70  mm  specimen  data  shown  in  Figure  21  is  small,  it  is  likely  that  a  thicker  specimen  would 
not  produce  higher  crack  growth  rates  than  the  12.70  mm  thick  specimen,  indicating  that 
thickness-independent  data  is  probably  provided  by  the  12.70  mm  thick  specimen. 

Figure  22  shows  the  comparison  of  the  crack  growth  rates  for  6.35  mm  and  12.70  mm 
thick  specimens.  Both  tests  were  conducted  at  760°C  (1400°F),  120  sec  dwell,  R  =  0.1.  The 
crack  growth  data  from  the  thin  specimen  is  nearly  equal  to  the  thick  specimen  data. 
Therefore,  a  12.70  mm  thick  specimen  should  be  thick  enough  to  give  thickness-independent 
data. 


To  confirm  the  conclusion  that  12.70  mm  thick  specimens  produce  thickness-independent 
data  for  dwell  tests,  a  19.05  mm  thick  specimen  was  tested  at  760°C,  R  =  0.1,  600  sec  dwell  at 
maximum  tensile  load.  The  comparison  of  the  data  with  data  from  a  12.70  mm  thick  specimen 
tested  at  the  same  conditions  shown  in  Figure  23  indicates  that  12.70  mm  thick  specimens 
produce  thickness-independent  data. 

Note  that  the  thickness  required  at  649°C  to  assure  plane  strain  constraint  at  the  crack 
tip  increases  from  6.35  mm  for  cyclic  tests  to  12.70  mm  (or  possibly  greater)  for  tests 
containing  dwells  at  maximum  tensile  load.  One  explanation  of  this  behavior  can  be  obtained 
by  considering  the  environmental  influences  (oxidation)  on  crack  tip  elasticity.  Environmental 
degradation  promotes  crack  tip  embrittlement,  which  in  turn  acts  to  discourage 
through-thickness  strain.  Therefore,  thinner  sections  will  exhibit  more  enhanced  plane  strain 
conditions  in  an  oxidizing  atmosphere  than  in  an  unoxidized  condition.  During  sustained  load, 
or  long  dwell,  tests,  a  protective  oxide  layer  can  form  which  prevents  further  oxygen  diffusion 
into  the  material.  As  this  prophylactic  effect  reduces  the  oxidation  rate,  the  thickness  required 
for  plane  strain  constraint  at  the  crack  tip  increases. 

Constant  load  amplitude  fatigue  testing  to  determine  the  effects  of  thickness  on  crack 
propagation  indicates  that  all  cyclic  tests  should  be  performed  on  6.35  mm  (0.250  in.)  thick 
specimens,  and  cyclic/dwell  (and  sustained  load)  tests  should  be  performed  on  12.70  mm 
(0.500  in.)  thick  specimens  for  thickness-independent  data  at  the  temperatures  tested  (427°C, 
649°C,  760°C).  These  represent  minimum  thicknesses. 

Fraclography 

The  fracture  surfaces  of  specimens  1102  (427°C,  cyclic),  1101  (649°C,  cyclic),  and  1106 
(649°C,  cyclic)  were  examined  using  the  transmission  electron  microscope  (TEM).  Five  areas 
of  each  fracture  surface  were  scrutinized  to  determine  the  fracture  mechanism  at  different 
stress  intensity  (K)  levels  as  shown  in  Figures  24,  25  and  26.  A  more  in-depth  study  coupled 
with  a  scanning  electron  microscope  evaluation  might  be  more  conclusive,  but  is  not  within  the 
scope  of  this  program.  The  following  conclusions  were  determined  from  inspection  of  repre¬ 
sentative  areas. 


29 


All  specimens  had  a  heavy  oxide  on  the  fracture  surface,  which  obscured  most  fracture 
features.  Specimen  1101  showed  faint  oxidized  river  markings  in  area  1  (Figure  27),  but  the  other 
areas  were  too  oxidized  to  see  any  significant  features.  In  general  the  fracture  face  was  flat  and 
showed  no  signs  of  intergranular  propagation. 

There  was  also  evidence  of  oxidized  river  markings  in  areas  1  and  2  of  specimen  1106 
(Figure  28).  Area  3  (Figure  29)  appeared  to  have  remnant  striations,  but  the  heavy  oxide 
prevented  definite  identification.  This  fracture  surface  also  was  flat  and  showed  no  evidence  of 
intergranular  propagation. 

Specimen  1102  (427°C)  showed  definite  fatigue  striations  in  area  3  (Figure  30)  with  an 
average  spacing  of  6.9  X  10  4  mm  (AK  =  36)  which  correlates  well  with  the  da/dN  plot.  Area  4 
was  suggestive  of  remnant  coarse  striations  but  no  real  conclusion  could  be  made  from  what  little 
was  observed.  Figures  31  and  32  show  that  areas  1  and  2  were  relatively  flat  and  oxidized,  and  like 
the  other  two  specimens,  showed  no  evidence  of  intergranular  propagation. 
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Figure  24.  Macrofractograph  of  Specimen  1102  Showing  the  Areas 
Investigated 
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Figure  25.  Macrofractograph  of  Specimen  1101  Showing  the  Areas 
Investigated 
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Figure  26.  Macrofractograph  of  Specimen  1106  Showing  the  Areas  In¬ 
vestigated 
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Figure  28.  Compact  Specimen  No.  1106 
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Figure  29.  Compact  Specimen  So.  1106  Area  3 
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Figure  30.  Compact  Specimen  No.  1102  Area  3 
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Phase  II:  Development  ot  An  Interpolate  Model 

Introduction 

The  objective  of  Phase  II  was  to  develop  an  interpolate  model  to  permit  accurate 
predictions  of  crack  propagation  rates  for  an  advanced  powder  metallurgy  (P/M)  superalloy 
(AF2-1DA)  under  different  stress-temperature-time  conditions.  This  generalized  interpolate 
model,  based  on  the  hyperbolic  sine  equation 

log  (da/dN)  =  C,  sinh  (C2  (log  (AK)  +  C3))  +  C«,  (1) 

completely  describes  the  effects  of  stress  ratio  (R),  cyclic  frequency  (v  ),  and  tempe  ature  (T) 
on  the  crack  growth  rate  in  this  material. 

Interpolate  capacity  forms  the  fundamental  strength  of  the  hyperbolic  sine  model.  The 
procedure  known  as  the  interpolation  algorithm  for  calculating  the  SINH  coefficients,  describ¬ 
ing  fatigue  crack  propagation  (FCP)  for  any  frequency,  stress  ratio,  and  temperature,  is 
illustrated  in  the  following  paragraphs  (reprinted  from  Reference  7). 

Determination  of  SINH  Coefficients 

The  coefficients  (e.g.,  C2,  C3,  and  C,)  at  any  intermediate  value  of  an  element 
life-controlling  parameter,  can  be  determined  from  equation  9: 

C,  =  C*_  +  AC,;  j  =  2,3,4  (9) 

where: 

r Cj  i 

Cj  =  C3  J  =  interpolated  values  of  coefficients 

and: 

r  ACi  i 

AC,  =  I  AC3  I  =  differences  from  baseline  values 

L  AC,  J 


Since  the  SINH  coefficients  are  linear  functions  of  the  controlling  parameters,  *  it  is 
evident  that: 


AC,  1 

r  dC,/dv,  dCJdR,  dCJdT 

1 

p  Av 

ac3 

1  dC Jdv,  aCj/dR,  5CS/3T 

1  X 

AR 

AC,  J 

“  Lac  Jdv,  ac,/aR,  ac,/aT  . 

J  x 

L  AT 

NX1 

NXN 

NX1 

where: 


] 


(10) 


Av  "| 

AR  I  =  Differences  from  baseline  values 
AT  J 


'Strictly  speaking,  the  coefficents  are  nonlinear  functions  of  r,  E,  and  T;  however,  they  are  linear  functions  of  other 
functions.  This  simplification  was  used  here  for  presentation  clarity. 
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and  the  N  X  N  partial  derivative  matrix  is  easily  determined  from  the  slopes  of  the  lines  relating 
each  coefficient  with  each  rate  controlling  parameter.  The  computation  of  the  intermediate 
coefficients,  using  equation  9,  is  then  straightforward. 

Advanced  Regreaaion  Conalderatlona 

The  output  of  the  interpolation  algorithm  is  the  set  of  SINH  coefficients  describing  Fatigue 
Crack  Propagation  (FCP)  under  the  conditions  input.  A  parallel  consideration,  necessary  in  the 
modeling  phase,  is  some  multiple  regression  capability  which  will  allow  simultaneous 
consideration  of  several  different  collections  of  data,  each  differing  from  the  others  bv  only  one 
FCP  controlling  parameter.  The  data  for  one  test  condition  are  usually  regressed  separately:  data 
at  one  condition  are  not  allowed  to  influence  the  model  at  another  condition.  However,  if  the  final 
model  is  to  be  truly  interpolative.  behavior  at  one  condition  will  be  used  to  describe  FCP  at 
another  condition.  Therefore,  the  data  must  be  permitted  to  exhibit  their  mutual  influence 
during  the  modeling  process. 

Pratt  &  Whitney  Aircraft  has  developed  a  mathematical  technique  to  accomplish  this. 
Individual  sets  of  data  are  treated  independently  relative  to  some  of  the  SINH  coefficients,  while 
the  entire  collection  is  treated  as  an  entity  with  respect  to  the  interpolative  coefficents  (functions 
of  v.  R  or  T).  This  improved  model  has  been  programmed  (MOD-1)  and  is  in  use.  The  P&WA 
modeling  philosophy,  the  basic  model  formulation  including  characteristics  of  the  sinh  (math 
function)  and  SINH  (model),  as  well  as  the  SINH  descriptions  of  basic  FCP.  have  been  discussed 
previously.  The  actual  computational  procedure  required  to  perform  the  desired  data  modeling, 
referred  to  as  the  Method  of  Least  Squares,  is  described  in  detail  here. 

The  goal  of  this  procedure  is  to  determine  model  coefficients  so  that  the  resulting  curve 
through  the  data  will  have  the  least  (minimum)  summed  squared  error  between  calculated  and 
observed  values  for  the  dependent  variable  (Figure  ,13).  In  this  instance,  the  independent  and 
dependent  variables,  x  and  v,  are  log  (AK)  and  log  (da/dN),  respectively. 

Define  the  sum  of  the  squared  errors  as 

n 

EJ  =  Z  F,1  =  I  (y,,,,.  -  v,)’  (ID 

U  1 

Since  y,.al|  =  f(C,,  C,.  C4,  x,h  E  is  also  a  function  of  C2.  Cs.  C4. 


Now.  F,2  will  be  a  minimum  when  each  of  its  partial  derivatives  is  zero  simultaneous!  v.  Thai 
is 
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when  f  is  the  SINH  model. 


(ir.i 

(Hit 

(17) 


Now,  substituting  equations  15,  16.  17  and  18  into  equations  12.  l.'l  and  14.  and  solving*  the 
resulting  three  simultaneous  equations  provides  the  values  for  C2.  C3  and  C,  for  which  equation 
11  will  be  a  minimum. 
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Figure  33.  Method  of  Least  Squares 


*  In  this  instance,  the  resulting  simultaneous  equations  are  nonlinear  in  and  C4.  The  solution  therefore  requires 

some  iterative  procedure  such  as  an  N -dimensional  Newton -Raphson  method. 


The  foregoing  discussion  explains  the  procedure  for  determining  the  coefficients  for  one 
SINH  curve.  Suppose  further  that  each  SINH  representation  of  FCP  is  related  to  each  other, 
the  relationship  depending  on  differences  in  frequency,  stress  ratio,  or  temperature.  Consider 
frequency  as  an  example,  and  assume  that  the  points  of  inflection  are  linearly  related,  viz: 

C:,J  =  C„  f  C*  (C„j)  (19) 

for  j  different  SINH  curves. 

Assume  also  that  C2  and  C,  are  related  to  test  frequency  by  equations  20  and  21. 

c„j  =  +  cM  log  (.-)  (20) 

C„j  =  C37  +  CM  log  (e)  (21) 

Coefficients  CM  through  CM  can  then  be  determined  by  substituting  equations  19,  20  and  21 
into  equation  15  and  differentiating  with  respect  to  C:„  through  Cw  in  a  manner  analogous  to 
that  used  in  determining  C2,  C3,  and  C,  in  the  foregoing  discussion. 

Effect  of  Cyclic  Frequency 

Previous  experience  with  P/M  alloys  indicates  that  changing  test  frequency,  while 
holding  stress  ratio  and  temperature  constant,  produces  crack  growth  curves  similar  in  shape 
but  shifted  along  a  nearly  vertical  line  passing  through  the  points  of  inflection.  The  location  of 
these  inflection  points  is  related  to  test  frequency:  reduced  frequency  (increased  cycle 
duration)  are  observed  to  accelerate  crack  growth  rate. 

The  AF2-1DA  crack  growth  rates  for  cyclic  frequencies  0.0082  Hz,  0.17  Hz,  and  20  Hz 
(649°C ,  R  =  0.5)  are  shown  individually  in  Figures  24,  25  and  26,  respectively.  The 
comparison  of  these  crack  growth  rates  is  shown  in  Figure  27;  note  that  points  of  inflection  are 
represented  by  a  piecewise  linear  function,  with  each  segment  having  the  form:  C,  =  C\,+  CM 
^  ^  4*  figure  28  illustrates  the  relationships  between  C„  ( 1 ,,  C4,  and  frequency. 

The  coefficients  are  linearly  related  and  equations  for  all  interpolative  SINH  frequency 
model  coefficients  are  shown  in  Figure  28,  except  C,  which  was  held  constant  at  0.5.  This 
model  completely  describes  the  effect  of  frequency  on  FCP  in  AF2-1DA  at  this  temperature 
and  stress  ratio. 

Effect  of  Street  Ratio 

As  with  frequency,  stress  ratio  effects  are  easily  described  using  the  SINH  model.  The 
individual  data  sets  used  in  the  649°C,  20  Hz  stress  ratio  model  (R  =  0.1,  R  =  0.5,  R  =  0.8) 
are  shown  in  Figures  39,  40  and  41.  The  comparison  of  these  crack  growth  rates  is  shown  in 
Figure  42.  As  indicated  by  Figure  42,  it  was  necessary  to  use  a  piece-wise  linear  function  to 
describe  the  relationship  between  C3  and  C,.  Kach  segment  of  this  function  is  of  the  form: 

c,  =  c»  +  c„  X  ct. 

The  SINH  crack  growth  rate  curve  can  be  calculated  for  any  positive  stress  ratio  by  using  the 
relationships  between  C„  C„  and  C4,  and  stress  ratio  given  in  Figure  43.  Again,  C,  was  held 
constant  at  0.5  completing  the  model.  This  model  completely  describes  the  effects  of  stress 
ratio  at  this  temperature  and  frequency. 
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Figure  35.  AF2-1DA  Crack  Growth  Data  at  649°C,  0.17  Hz,  R  =  0.1 


AK  -  MPav/m 


AK  -  ksi\/  in. 


Figure  36.  AF2-1DA  Crack  Growth  Data  at  649°C,  20  Hz,  R  =  0.1 
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Figure  39.  AF2-IDA  Crack  Growth  Data  at  649°C,  20  Hz,  R  =  0.1 
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Figure  40.  AF2-1DA  Crack  Growth  Data  at  649cCt  20  Hz t  R  =  0.5 
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The  individual  data  sets  used  in  the  649°C,  0.17  Hz  stress  ratio  model  (R  =  0.1,  R  =  0.5) 
appear  in  Figures  44  and  45.  The  points  of  inflection  describe  a  straight  line  in  Figure  46  such 
that  C3  (horizontal  inflection  location)  is  a  linear  function  of  C4  (vertical  location);  C,  =  CM  + 
C34  X  C4.  The  SINH  crack  growth  rate  curve  can  be  calculated  for  any  positive  stress  ratio 
(keeping  temperature  and  frequency  constant)  by  using  the  relationship  between  stress  ratio 
and  C2,  C3,  and  C4  given  in  Figure  47. 

The  individual  data  sets  used  in  the  649°C,  0.0083  Hz  stress  ratio  model  (R  =  0.1, 
R  =  0.5,  R  =  0.8)  appear  in  figures  48,  49,  and  50.  Again,  comparison  of  these  crack  growth 
rates  in  Figure  51  indicates  that  the  points  of  inflection  describe  a  straight  line.  The  SINH 
crack  growth  rate  curve  can  be  calculated  for  any  positive  stress  ratio  by  using  the  relationship 
between  stress  ratio  and  C2,  C3,  and  C4  given  in  Figure  52. 

Since  the  effect  of  stress  ratio  on  crack  growth  rates  can  vary  significantly  with  frequency 
(Reference  1),  a  stress  ratio  model  is  needed  near  the  extremes  of  the  frequencies  of  interest. 
In  this  case,  the  20  Hz  model  and  the  0.0083  Hz  model,  which  were  discussed  in  the  preceding 
paragraphs,  are  used  in  conjunction  with  the  0.17  Hz  model  to  describe  the  effects  of  stress 
ratio  over  the  active  frequency  range. 

Effect  of  Temperature 

Similar  to  frequency  and  stress  ratio,  the  effects  of  temperature  can  be  described  using 
the  SINH  model.  Figures  53,  54  and  55  show  the  individual  data  sets  (427,  649  and  760°C) 
used  in  the  temperature  model.  The  composite  plot  of  these  data  sets  (Figure  56)  shows  that 
once  again  the  inflection  points  describe  the  straight  line  C,  =  Ca  +  Ca  X  C4.  The 
relationship  between  C2,  C4,  and  temperature  is  shown  in  Figure  57  with  the  equations  for  all 
the  SINH  model  coefficients  except  C,  (C,  =  0.5).  This  model  completely  describes  the  effects 
of  temperature  on  AF2-1DA  crack  growth  rates  even  where  data  are  unavailable. 

Effect  of  Hold  Time 

The  length  of  time  a  crack  tip  is  subjected  to  maximum  tensile  load  each  cycle  will 
obviously  affect  the  crack  growth  rate.  The  differences  in  propagation  rates  due  to  varying 
hold  times  can  be  characterized  by  the  SINH  model.  Figures  58  through  61  illustrate  the 
649°C,  R  =  0.1,  individual  data  sets  at  30  sec,  120  sec,  300  sec,  and  600  sec  hold  times, 
respectively.  The  effect  of  these  hold  times  on  crack  growth  rate  can  easily  be  seen  in  Figure 
62,  as  can  the  linear  relationship  (C3  =  +  CM  X  C4)  between  the  inflection  points.  The 

relationship  between  C2,  C4  and  hold  time  is  shown  in  Figure  63.  This  SINH  model  describes 
the  effects  of  maximum  tensile  hold  time  on  AF2-1DA  crack  growth  rates  between  30  sec  and 
600  sec  hold  times  (649°C,  R  =  0.1),  even  where  data  are  unavailable. 

As  temperature  increases,  the  amount  of  oxidation  and  creep  at  the  crack  tip  also 
increases,  therefore,  larger  differences  in  crack  growth  rates  due  to  changes  in  hold  times  are 
expected.  The  individual  data  sets  for  the  7606C  dwell  model  are  shown  in  Figures  64  and  65 
(zero  and  120  sec  dwells).  Only  a  two-point  model  is  used  (Figure  66)  for  reasons  explained  in 
the  Model  Verification  Section.  The  SINH  crack  growth  rate  curve  can  be  calculated  for  any 
hold  time  between  zero  and  120  seconds  (760°C,  R  =  0.1)  by  using  the  relationships,  given  in 
Figure  67,  between  C2,  C„  C4  and  maximum  tensile  hold  time. 
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Figure  44.  AF2-1DA  Crack  Growth  Data  at  649°C.  0.17  Hz,  R  =  0.1 
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Figure  45.  AF2-1DA  Crack  Growth  Data  at  649°C,  0.17  Hz,  R  =  0.5 
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Coefficients  C2  and  C4  vs  (1  -  Stress  Ratio) 
C2  =  6.2319  +  1.6606  Log  (1  -  R  Ratio) 

C4  -  -4.3939  +  1.6845  Log  (1  -  R  Ratio) 
C3  =  -4.3173  -  0.6326  C4 
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Figure  17.  Effect  of  Sires : 

s-  Ratio  on  SINH  Coefficients,  649°C , 
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Figure  49.  AF2-IDA  Crack  Growth  Data  at  649°C,  0.0083  Hz,  R=  0.5 
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Figure  50.  AF2-IDA  Crack  Growth  Data  at  649°C.  0.0083  Hz,  R—  0.8 
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Figure  52.  Effect  of  Stress  Ratio  on  SINH  Coefficients,  649°C,  0.0083  Hz 
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Figure  53.  AF2-1DA  Crack  Growth  Data  at  427eC,  0.17  Hz,  R  =  0.1 
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gure  54.  AF2-1DA  Crack  Growth  Data  at  649°C ,  0.17  Hz,  R  «  0.1 
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Figure  55.  AF2-1DA  Crack  Growth  Data  at  760°C,  0.17  Hz,  R  «  0.1 
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Figure  60.  AF2-1DA  Crack  Growth  Data  at  649°C,  300  Sec  Dwell,  R  =  0.1 
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Figure  61.  AF2-1DA  Crack  Growth  Data  at  649°C,  600  Sec  Dwell,  R  -  0.1 
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Fieure  62.  Effect  of  Maximum  Tensile  Hold  Time  on  AF2-1DA  Crack  Growth  Rate  at 
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Coefficients  C2  and  C4  vs  Dwell  Time 
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Figure  63.  Effect  of  Dwell  Time  on  SINH  Coefficients  at  649° C,  R  =  0.1 
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Phase  III:  Model  Demonstration 


Statistical  Confidence  of  the  Model 

A  mathematical  model  describing  crack  propagation  has  utilitarian  value  only  insofar  as  it 
can  be  used  in  life  prediction;  the  overall  accuracy  of  a  model  can  be  measured  by  the  accuracy 
of  the  resulting  prediction.  To  provide  a  basis  for  comparing  the  accuracy  of  various  life 
predictions,  a  simple  correlative  parameter  is  used,  N1,r,.,1/N„l.t,  the  quotient  of  predicted  and 
actual  cyclic  lives.  Ideally,  this  quotient  is  1.0  and  decimal  deviations  from  the  ideal  can  be 
quickly  interpreted  as  percent  error  of  the  prediction. 

One  of  the  objectives  of  model  development  is  to  quantify  the  intrinsic  error  of  the  model. 
The  fundamental  test  of  model  validity  is  its  ability  to  “predict"  the  lives  of  the  specimens  which 
were  actually  used  in  its  development.  Therefore,  the  ex  post  facto  lives  of  all  specimens 
contributing  to  current  model  development  have  been  calculated  and  statistics  of  model  success 
collected.  The  resulting  distribution  of  N„r,.,i/Nm,  provides  a  quantitative  assessment  of  the 
prediction  error  to  be  expected  during  the  model  demonstration. 

The  procedure  used  for  characterization  of  the  statistical  reliability  of  the  crack  propagation 
models  is  summarized  as  follows: 

•  Generate  an  interpolative  crack  propagation  model  from  a  series  of  related 
(da/dN,  AK)  data  sets  using  multiple  regression  capability  (MOD-1). 

•  Using  the  generated  model,  predict  the  (da/dN.  AK)  curves  for  the  test 
conditions  of  the  specimens  used  to  generate  the  model. 

•  Produce  specimen  life  calculations  by  integrating  the  crack  propagation 
curves. 

•  Calculate  a  value  of  N„r„,/Nm.,  for  each  of  the  test  specimens  employed  in  the 
model. 

•  After  completing  the  above  four  steps  for  each  separate  model  (e.g..  cyclic 
frequency,  stress  ratio,  temperature),  the  entire  population  of  values  of 
N1)r,.<t/Nm.,  for  all  models  is  analyzed  statistically. 

•  The  results  of  the  statistical  analysis  are  used  to  characterize  model 
reliability  and  suggest  possible  model  improvements. 

This  process  has  been  applied  to  all  models  presented  in  this  report,  and  the  resulting 
statistics  have  been  a  significant  aid  to  model  evaluation.  Such  an  assessment  fully  describes  the 
error  in  the  model,  since  linear  relationships  between  similar  SINH  coefficients  are  imposed.  By 
requiring  that  the  linear  model  functions  precisely  locate  the  modeling  curves,  a  rigid 
interpolation  model  is  obtained,  and  all  error  in  the  model  is  forced  to  appear  in  the  SINH  "fits" 
of  the  individual  data  sets.  Therefore,  an  analysis  of  the  collection  of  values  of  Npr.VN,,,., 
associated  with  a  particular  model  is  a  critical  test  of  the  inherent  error  of  the  model. 

The  results  of  the  statistical  analysis  of  each  of  the  crack  propagation  models  presented 
in  this  report  are  found  in  Table  4.  The  statistical  parameters  derived  from  each  model 
represent  calculations  based  on  an  assumed  log-normal  probability  distribution  for  all  such 
data.  The  collection  of  values  of  Npre<(/Na(.t  for  all  specimens  used  to  generate  the  models  is 
presented  in  a  log-normal  probability  plot  in  Figure  68.  The  mean  value  of  this  sample 
population  is  0.975  and  the  log  standard  deviation  is  0.0514.  This  implies  that  the  error  in  life 
prediction  may  be  expected  to  be  less  than  25  '7  for  95 cc  of  the  future  life  predictions  using 
the  models. 
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Experimental  Verification 


Two  compact  specimens  were  tested  at  each  of  two  previously  untried  conditions  to  verify 
the  interpolative,  predictive  ability  of  the  SINH  crack  growth  model.  The  two  test  conditions 
(538°C,  R  =  0.65,  0.017  Hz  and  718°C,  R  =  0.1,  15-sec  dwell  at  maximum  tensile  load)  were 
provided  by  the  Air  Force  Project  Engineer.  The  predicted  crack  growth  rate  curves  were  given 
to  him  prior  to  testing  with  an  explanation  for  revising  the  dwell  model  to  predict  the  15-sec 
dwell  test. 

Revision  of  the  dwell  model  was  required  since  prediction  of  the  verification  test  crack 
growth  rate  required  extrapolation  into  a  region  where  crack  extension  is  not  a  linear  function 
of  dwell  time.  Figure  69  depicts  the  effects  of  dwell  time  on  crack  growth  rates  at  760°C.  The 
original  interpolative  dwell  model  employed  a  minimum  dwell  time  of  120  sec.  However,  one  of 
the  verification  tests  included  a  15  sec  dwell.  As  indicated  by  Figure  69,  a  linear  extrapolation 
from  120  sec  dwell  to  15  sec  dwell  (dotted  lines)  would  result  in  a  significant  deviation  from 
the  data.  Therefore,  the  dwell  model  was  revised  to  take  0  to  120  sec  dwell  times  into  account. 
The  revision  allows  the  model  to  be  interpolated  between  0  and  120  sec  dwell  times,  rather 
than  extrapolating  from  the  120  sec  dwell,  resulting  in  a  more  accurate  representation  of  actual 
conditions.  All  dwell  predictions  were  made  using  the  revised  model. 


Figure  70  presents  a  comparison  of  the  da/dN  vs  AK  curve  and  the  actual  data  for  the 
538°C,  R  =  0.65,  0.017  Hz  verification  testing.  The  curves  correlate  well  with  the  data, 
Npred/N.ct  =  0.789  for  Specimen  1167  and  Npred/Nact  =  1.008  for  Specimen  1168.  Therefore, 
the  specimen  life  predictions  detailed  in  Figures  71  and  72  fall  well  within  the  expected  limits 
of  the  model  (i.e.,  less  than  25%  error). 

Figure  73  illustrates  the  comparison  of  the  predicted  da/dN  vs  AK  curve  and  the  actual 
data  for  the  718°C.  R  =  0.1,  15  sec  dwell  testing.  Here  again,  the  curves  describe  the  data  well, 
Npred/Nact  =  0  843  for  Specimen  1153  and  Npred/Nact  =  1.022  for  Specimen  1154.  Thus, 
specimen  life  predictions  shown  in  Figures  74  and  75  fall  within  the  expected  limits  of  the 
model  previously  described. 

Conciuaiona 

•  The  interpolative  hyperbolic  sine  (SINH)  model  accurately  predicts  crack 
propagation  life  of  AF2-1DA  at  elevated  temperatures. 

•  At  elevated  temperatures  (427°C  <  T  <  760°C),  the  thick¬ 

ness-independent  data  require  specimens  for  cyclic  tests  to  be  a  minimum 
of  6.35  mm  thick  and  for  cyclic/dwell  tests  to  be  a  minimum  of  12.70  mm 
thick. 

•  Additional  studies  are  needed  at  various  temperatures  to  determine  the 
effect  of  dwell  length  on  crack  growth  rates. 


TABLE  4. 

STATISTICAL  EVALUATION  OF  THE  CRACK 
PROPAGATION  MODELS 


20  Hz,  1200° F  Stress  Ratio  Model 


Spec  No. 

1125 

0.1 

0.930 

1135 

0.5 

0.982 

1142 

0.8 

1.036 

Mean  Np/NA 

-  0.982 

Log  (»)  - 

0.0083  Hz, 

1200° F  Stress  Ratio  Model 

Spec  No. 

R^ 

Mu 

1104 

0.1 

1.031 

1115 

0.5 

1.067 

1144 

0.8 

1.436 

Mean  Np/NA 

-  1.164 

Log  (e)  - 

Frequency  Model 

Spec  No. 

njna 

1104 

0.1 

0.923 

1101 

0.1 

0.930 

1126 

0.1 

1.002 

1125 

0.1 

1.031 

Mean  Np/NA 

-  0.970 

Log  (o)  - 

I400°F  Dwell  Model 

Spec  No. 

R_ 

njna 

1149 

0.1 

0.825 

1148 

0.1 

1.027 

Mean  N,/NA 

-  0.920 

I<og  (o)  = 

10  cpm,  h 

!  —  0.1  Temperature  Model 

Spec  No. 

JL 
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0.1 

0.999 

1120 

0.1 

0.915 

1101 

0.1 

1.002 

1126 

0.1 

0.923 
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0.1 

0.825 

Mean  Np/NA 

-  0.930 

Log  (<r)  - 

1200°F  Dwell  Model 

Spec  No. 

R^ 

n„/na 

1165 

0.05 

0.761 

1103 

0.1 

0.922 

1152 

0.1 

1.037 

1124 

0.1 

1.042 

Mean  Np/NA 

-  0.930 

Log  (it)  - 

10  cpm, 

1200°  F  Stress  Ratio  Model 

Spec  No. 

R^ 

Nn'NA 
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0.1 
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0.1 
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Figure  71.  Life  Prediction  for  Verification  Specimen  1167  Tested  at  .5.18° C. 

R  =  0.65.  0.017  Hz 
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Figure  72.  Life  Prediction  for  Verification  Specimen  1168  Tested  at  538°C, 

R  =  0.65,  0.017  Hz 
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Figure  74.  Life  Prediction  for  Verification  Specimen  1 153  Tested  at  7 1 8°  C,  R  - 
O.I,  15  sec  Dwell 
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Figure  75.  Life  Prediction  for  Verification  Specimen  1154  Tested  at  718°C,  R  = 
0.1,  15  sec  Dwell 
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